Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 



CP Violation and Nonleptonic B -meson Decays 

Y.-Y. Keum and A. I. Sanda 

Department of Physics Nagoya University 



o 
o 

(N 



We discuss the perturbative QCD approach for exclusive non-leptonic two body B-meson decays. We briefly review its ingredients 
and some important theoretical issues related to the factorization approach. PQCD results are compatible with present experimental 
data for the charmless B-meson decays. We predict the possibility of large direct CP violation effects in B° — » tz + k~ (23 ±7%) and 
B° — » K + n~ ( — 17 ±5%). For charmed decay B — > D^K, we get large non-factorizable contributions for Color-suppressed decays and 
obtain \d2/a\ | ~ 0.4 — 0.5 and Argia^jai ) ~ —42°. In the last section we investigate the method to extract the weak phases fa from 
B — » nn decay and fa from Kit modes. From BaBar measurement of CP asymmetry for the K + n~ decay, the prefered CKM weak 
phases are: fa = (24 ±2)°, fa = (78 ±22)° and fa = (78 ±22)°. 
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1 Introduction 

Understanding nonleptonic B meson decays is crucial for 
testing the standard model, and also for uncovering the 
trace of new physics. The simplest case is two-body non- 
leptonic B meson decays, for which Bauer, Stech and 
Wirbel proposed the factorization assumption (FA) in their 
pioneering work [[Q. Considerable progress, including the 
generalized FA [ 13] 0] and QCD-improved FA (QCDF) 
[ 0, has been made since this proposal. On the other 
hand, techniques to analyze hard exclusive hadronic scat- 
tering has been developed by Brodsky and Lepage [ 6| 
based on collinear factorization theorem in perturbative 
QCD (PQCD). A modified framework based on kj factor- 
ization theorem has been given in [[710> and extended to 
exclusive B meson decays in rl9l ll0l[TTIll2l . The infrared 
finiteness and gauge invariance of kj factorization theorem 
was shown explicitly in f 1 1 31 . Using the PQCD approach, 
we have investigated the dynamics of nonleptonic B meson 
decays r il4N15llT6l . Our observations are summarized as 
follows: 

1. FA is approximately correct, as our computation 
shows that nonfactorizable contributions in charm- 
less B meson decays are negligible. 

2. Penguin amplitudes are enhanced, as the PQCD for- 
malism inludes dynamics from the energy region as 
low as \Jhm,b. Here A = mj — nih, being the B me- 
son and b quark mass difference. Note that yA% 
is much lower than mi,{or ra/,/2) which is often used 
as the energy scale. 

3. Annihilation diagrams give rise to large short- 
distance strong phases through (S + P)(S — P) pen- 
guin operators. 



4. The sign and magnitude of CP asymmetries in two- 
body nonleptonic B meson decays can be calculated, 
and we have predicted relatively large CP asymme- 
tries in the B — ► K^X [[T4| and nit modesr lT5llT6l . 

In this talk we summarize the PQCD method and important 
theoretical issues, and describe the analysis of branching 
ratios of B-meson decays including B — > D^*'7C decays and 
discuss the origin of large direct CP-violation in B — > 
process. 

2 Ingredients of PQCD and Theoretical Is- 
sues 

End Point Singularity and Form Factors: If we calculate 
the B — > % transition form factor F Blz at large recoil using 
the Brodsky-Lepage formalism M171 1181 . a difficulty im- 
mediately occurs. The lowest-order diagram for the hard 
amplitude is proportional to \j{x\X^), x\ being the mo- 
mentum fraction associated with the spectator quark on the 
B meson side and xt, the momentum fraction carried by 
the spectator quark in final pion. If the pion distribution 
amplitude vanishes like X3 as X3 — > (in the leading-twist, 
i.e., twist-2 case), F Bk is logarithmically divergent. If the 
pion distribution amplitude is a constant as X3 — > (in the 
next-to-leading-twist, i.e., twist-3 case), F Bk even becomes 
linearly divergent. These end-point singularities have also 
appeared in the evaluation of the nonfactorizable and anni- 
hilation amplitudes in QCDF. 

Note that in the above discussion parton transverse mo- 
menta k± has been neglected. More accurately, we have 
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(x 3 M%+q x ) [ Xl x 3 Mi+(k^-k 3± y 



(i) 



and the end-point singularity is smeared out. More precise 
analysis including the Sudakov and threshold resummation 
effects has been given [ 14 1. 
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In PQCD, we can calculate the form factors for B — ► P, V 
transitions r fT9ll20l . 

Dynamical Penguin Enhancement vs Chiral Enhance- 
ment: The typical hard scale is about 1.5 GeV as discussed 
in Ref. N14l . Since the RG evolution of the Wilson coeffi- 
cients Cnfi(t) increase drastically as t < Mb/2, while that 
°f C\ %{t) remain almost constant, we can get a large en- 
hancement effects from both Wilson coefficents and matrix 
elements in PQCD. 

In general the amplitude can be expressed as 

Amp ~ [«i,2 ± fl4 ± m^ V (jJ,)a6] ■ < Kn\0\B > (2) 

with the chiral factors m^^i) — mp/[mi(jx) +1112(11)] for 
pseudoscalar meson and mi = my for vector meson. To ac- 
commodate the B — » Kn data in FA and QCD-factorization 
approach, one relies on the chiral enhancement by increas- 
ing mo to as large as about 3 GeV at /I = m& scale. So two 
methods accomodate large branching ratios of B — > Kn. 
But, there is no such adjustable parameter for B — > PV de- 
cay. For B — > PV there is no chiral factor in LCDAs of the 
vector meson. Here mo = my. It is difficult to explain large 
penguin contribution in QCDF[|22). 

We can test whether a dynamical enhancement or a chiral 
enhancement is responsible for the large B — > Kn branch- 
ing ratios by measuring the B — ► VP,VV modes. In these 
modes penguin contributions dominate, such that their 
branching ratios are insensitive to the variation of the uni- 
tarity angle 03 . Our prediction for various modes are shown 
at Table 2. We point out that QCDF can not globally fit the 
experimental data for B — > PP,VP and VV modes simul- 
taneously with same sets of free parameters (ph,<Ph) ar >d 

Fat Imaginary Penguin in Annihilation: There is a folk- 
lore that annihilation contribution is negligible compared 
to the W-emission. For this reason, annihilation contribu- 
tion was not included in the general factorization approach 
and the first paper on QCD-factorization by Beneke et al. 
r 1241 . In fact there is a suppression effect for the opera- 
tors with structure (V — A) (V — A) because of a mechanism 
similar to the helicity suppression for n —> jJ-V^. How- 
ever annihilation from the operators 05,6,7,8 with the struc- 
ture (S — P)(S + P) via Fiertz transformation survive un- 
der the helicity suppression. Moreover, they provide large 
imaginary part. The real part of the factorized annihilation 
contribution becomes small because there is a cancellation 
between left-handed gluon exchanged diagram and right- 
handed gluon exchanged one as shown in Table 1 of ref.[ 
1161 . This mostly pure imaginary value of annihilation is a 
main source of large CP asymmetry in n + n~ and K + n~ 
decays. In Table 3 we summarize the CP asymmetry in 
B — > K(n)n decays with experimental measurements. 



Small Strong Phase for FA and QCDF: We have seen 
that the dominant strong phase in PQCD comes from the 
factorizable annihilation diagram[ 14j . For FA and QCDF, 
stong phases come from the Bander-Silverman-Soni (BSS) 
mechanism[ 21 1 and from the final state interaction (FSI). 
In fact, the two sources of strong phases in the FA and 
QCDF are strongly suppressed by the charm mass thresh- 
old and by the end-point behavior of meson wave func- 
tions. So the strong phase in these approaches is almost 
zero without soft-annihilation contributions. 



3 Numerical Results in Charmless B-decays 

Branching ratios in B -> PP,VP and VV: The PQCD ap- 
proach allows us to calculate the amplitudes for charmless 
B-meson decays in terms of light-cone distribution ampli- 
tudes upto twist-3. We focus on decays whose branch- 
ing ratios have already been measured. We take allowed 
ranges of shape parameter for the B-meson wave function 
as cob = 0.36 — 0.44, which accommodate reasonable form 
factors: both F Bn {0) = 0.27 - 0.33 and F BK (0) = 0.31 
0.40. We use values of chiral factor with mfi = \3GeV 
and Wq = 1 .IGeV . Finally we obtain branching ratios for 
B -> K(n)n [ Ql E3, K<j> [ [19] |25) K*<j>[ |26| and K*n[ 
1141 . which are in agreement with experimental data. In 
table 1 and 2 we summarize our predictions for various de- 
cays with updated experimental measurements which was 
reported at FPCP-2003 workshop M27l. 

CP Asymmetry of B — > nn,Kn: Because we have a 
large imaginary contribution from factorized annihilation 
diagrams, we predict large CP asymmetry (~ 25%) in 
B° — > n + n~ decay and about —15% CP violation effects 
in B° — » K + n~. The detail prediction is given in Table 
3. The precise measurement of direct CP asymmetry (both 
magnitude and sign) is a crucial way to test factorization 
models which have different sources of strong phases. 

Understanding of Br(K*n) and Br(coK + / con + ): In 
PQCD, penguin contributions enhances both Br(Kn) and 
Br(K*n). As noted before, in FA and QCDF, the penguin 
enhancement can achieved by taking niQ as large as 3 GeV , 
however, they fail to explain the large branching ratio for 
K*n decay. 

Another hot issue is how can we understand 
Br(coK/con) ~ 1. Since coK + decay is penguin dom- 
inant process and con + decay is tree dominant one: 
(a\ + xa2), it is hard to get the same branching ratio. 
PQCD method predicts Br(a>K + ) - 3.22 x 10~ 6 and 
Br(a>n + ) — 6.20 x 10~ 6 , which still has about factor 2 
differences between them. We need more precise measure- 
ments on coK + decay. 



Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 



3 




Figure 1. Plot of A„ K versus S x% for various values of <j) 2 with 0i = 24.3°, 0.18 < R c < 0.30 and -41° < 8 < -32° in the pQCD 
method. 



4 Extraction of 02 from B — > n + n 

Even though isospin analysis of B — > 7T7T can provide a 
clean way to determine 02, it might be difficult in practice 
because of the small branching ratio of B° — > 7t°n°. Here 
we describe the time-dependent analysis of B°(t) — > 7T + 7C~. 
Since penguin contributions are sizable and are about 20- 
30 % of the total amplitude, we expect that direct CP vio- 
lation can be large if strong phases are different in the tree 
and penguin diagrams. 

The ratio between penguin and tree amplitudes is R c = 
\P C JT C \ (here we use the c-convention notation) and the 
strong phase difference between penguin and tree ampli- 
tudes 8 = 8p — St- The time-dependent asymmetry mea- 
surement provides two equations for C nK and S KK in terms 
of three unknown variables R c , 8 and 02 [ 28 1. 

Since pQCD provides us Rc = and -41° < 8 < 

—32°, the allowed range of 02 at present stage is deter- 
mined as 55" < 02 < 100° as shown in FigurefJ] 

According to the power counting rule in the pQCD ap- 
proach, the factorizable annihilation contribution with 
large imaginary part becomes big and give a negative 
strong phase from —in8(k 2 L — xMg). Therefore we 
have a relatively large strong phase in contrast to QCD- 
factorization (8 ~ 0") and predict large direct CP violation 
effect in B Q -> n + 7I~ wifhA ep (fl° -> n + 7t~) = (23 ±7)%, 
which will be tested by more precise experimental mea- 
surement within two years. 



In the numerical analysis, since the data by Belle collabo- 
ration!^ 29 1 is located outside allowed physical regions, we 
considered the recent BaBar measurement! 30 1 with 90% 
C.L. interval taking into account the systematic errors: 

• S x% = 0.02 ±0.34 ±0.05 [-0.54, +0.58] 

• A K7l = 0.30 ±0.25 ±0.04 [-0.72, +0.12]. 

The central point of BaBar data corresponds to 0i = 24°, 

02 = 78°, and 3 = 78" when the Standard Model works. 

Even though the data by Belle collaboration [ 29 1 is located 
outside allowed physical regions, we can have overlapped 
ranges within 2 a bounds. 

5 Extraction of <fo from B — > Kk 

By using tree-penguin interference in B° — > K + n~(~ T + 
P 1 ) versus B+ -> K°7Z + {~ P'), CP-averaged B — > Kk 
branching fraction may lead to non-trivial constaints on the 

03 angle[ |3"T1 . In order to determine 03, we need one more 
useful information on CP-violating rate differencesT 1321 . 
Let's introduce the following observables : 



Rk 
A 



Br(B° — > K + n~) T 4 



Br(B + — > K°n^ 



= 1 — 2 tk cos8 CO503 + r| 



T(B° 



To 

-r(s° 



■K+7C- 



T(B- -►£%-) 



R K = 



^K°n+) 
-2rKsintyi sin8. 



(3) 
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Figure 2. Plot of R K versus <fe with r K = 0.164,0.201 and 0.238. 



where r% = \T' /p'\ is the ratio of tree to penguin ampli- 
tudes in B — > Kk and 8 = 8ji — 8pt is the strong phase dif- 
ference between tree and penguin ampli tides. After elimi- 
nating sin8 in Eq.(8)-(9), we have 



6 Large Nonfactorizable contribution in 
Charmed B-decays 



R K = 1 +r|± ^(4r K cos 2 fo-A 2 cot 2 fo). (4) 

Here we obtain r K = 0.201 ± 0.037 from the PQCD anal- 
ysis and Aq = —0.11 ±0.065 by combining recent BaBar 
measurement on CP asymmetry of B° — » K + n~: A cp (B° — ► 
K+n-) = -10.2±5.0±1.6% rl30l with present world av- 
eraged value offlr = 1.10±0.15rl33l. 

PQCD method provides 8 P > = 151", 8 T , = 1.4° and the 
negative cos8: cos8 = —0.91. As shown in Fig. 2, we can 
constrain the allowed range of 03 within 1 a range of World 
Averaged Rk as follows: 

• For cosS < 0, r% = 0. 164: we can exclude 0° < <fo < 
6°. 

• For cos8 <Q, rx — 0.201 : we can exclude 0° < 03 < 
6° and 35° < fa < 51°. 

• For cos8 < 0, rK = 0.238: we can exclude 0° < 03 < 
6° and 24° < 3 < 62°. 

When we take the central value of rK — 0.201, 03 is al- 
lowed within the ranges of 51" < 03 < 90°, which is con- 
sistent with the results by the model-independent CKM-fit 
in the (p, tj) plane. 



Large «2/ai in B — > D^n Decays: Being free from the 
end-point singularities, all topologies of decay amplitudes 
of charmed decays can be computed in the PQCD ap- 
proach, including the nonfactorizable color-suppressed di- 
agram. This amplitude can not be calculated in the QCDF 
approach based on collinear factorization theorem because 
of the existence of the singularities. However, we found 
that this amplitude is crucial for explaining the observed 
B — > Own branching ratios, since it is not suppressed by 
the Wilson coefficient (proportional to Ci/iV c ), and pro- 
vides a large strong phase required by the isospin relation. 
The tree annihilation amplitude, also contributing to the 
strong phase, is not important. As stated above, we have 
predicted large strong phases from the scalar-penguin an- 
nihilation amplitudes, which are required by the large CP 
asymmetries in two-body charmless decays. The success 
in predicting the strong phases from the nonfactorizable 
color-suppressed amplitudes for the two-body charmed de- 
cays further supports kj factorization theorem. The con- 
clusion is that the short-distance strong phase is sufficient 
to account for the B — > D^n data. A long-distance strong 
phase from final-state interaction may be small, though 
it should exist. Finally we obtained a2/a\ ~ 0.4 — 0.5 
and Arg(a2/ai) ~ —42° by including annihilation contri- 
butionsrl34l. 
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7 Summary and Outlook 

In this talk we have discussed ingredients of PQCD ap- 
proach and some important theoretical issues with nu- 
merical results. The PQCD factorization approach pro- 
vides a useful theoretical framework for a systematic anal- 
ysis on non-leptonic two-body B-meson decays. Specially 
pQCD predicted large direct CP asymmetries in B° — » 
n + n~ ,K + n~, K^ri^ and K^n decays, which will be 
a crucial touch stone to distinguish PQCD approach from 
others in future precise measurement. 

We discussed two methods to determine weak phases 
02 and 03 within the pQCD approach through 1) Time- 
dependent asymmetries in B° — > n + n~, 2) B — > Kit pro- 
cesses via penguin-tree interference. We can get interest- 
ing bounds on 02 an d 03 from present experimental mea- 
surements. From BaBar measurement of CP asymme- 
try in n + n~ decay the prefered CKM weak phases are: 
0! = (24 ±2)°, 02 = (78 ±22)° and 3 = (78 ±22)°. 
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Table 1. Branching ratios of B -> mt.Klt&nd KK decays with <p 3 = 80°, R b = \fp 2 + T} 2 = 0.38. Here we adopted = 1.3 GeV, 
mf = 1.7 GeV and 0.36 < ft) B < 0.44. Unit is 10~ 6 . 



Decay Channel 


CLEO 


BELLE 


BABAR 


World Av. 


PQCD 


n + n~ 


A c+1.4+0.5 


4.4 ±0.6 ±0.3 


4.7 ±0.6 ±0.2 


4.6 ±0.4 


5.93-10.99 


7t+K° 


4fi +1.8+0.6 


5.3±1.3±0.5 


5.5+o 9 ±0 - 6 


5.3 ±0.8 


2.72-4.79 


K°7t° 


<4.4 


<4.4 


<3.6 


< 3.6 


0.33-0.65 




18.U_ 21 _ 09 


18.5±1.0±0.7 


17.9±0.9±0.7 


18.2±0.8 


12.67-19.30 




1 R 0+3.7+2.1 
18.S_33_ 18 


22.0±1.9±1.1 


20.0±1.6±1.0 


20.6 ±1.4 


14.43 - 26.26 




n q+2.4+1.2 
iZ ' y -2.2-l.l 


12.8±1.4l};^ 


12.8l};o±1.0 


12.8 ±1.1 


7.87-14.21 


K°n° 


19 O+4.0+1.7 
1Z.S_33_ 14 


12.6 ±2.4 ±1.4 


10.4± 1.5 ± 1.8 


11.5 ± 1.7 


7.92-14.27 




<0.8 


<0.7 


<0.6 


<0.6 


0.06 


^+£° 


< 3.3 


<3.4 


<2.2 


<2.2 


1.4 




< 3.3 


<3.2 


< 1.6 


< 1.6 


1.4 



Table 2. Branching ratios of B -»• QK^and K*% decays with 3 = 80°, R b = \/p 2 + T] 2 = 0.38. Here we adopted = 1.3 GeV and 
mf = 1.7 GeV. Unit is lO" 6 . 



Decay Channel 


CLEO 


BELLE 


BABAR 


World Av. 


PQCD 


((>K 


5.5^ ±0.6 
5.4^ ±0.7 


9.4±1.1±0.7 
9.0±2.2±0.7 


10.0+^ ±0.5 
7.6ljj±0.5 


9.3±0.8 
7.7± 1.1 


8.1-14.1 
7.6-13.3 


<j)K*° 


10 6 +b - 4+L8 
iu.o_ 4 9 _j 6 

i i c+4.5+1.8 

1 1 ■ J -3.7-1.7 


r 7 +2.1+0.7 
-1.9-1.0 

io o +1 - 6+0 - 7 


12.1+^±1.1 
ll.l+Ji±0.8 


9.4±1.6 
10.7 ±1.1 


12.6-21.2 
11.5-19.8 




7-6^ ±1.6 


IO 4+4.2+4.1 
iV -^-3.9-7.1 


15.5±3.4±1.8 


12.3 ±2.6 


10.2-14.6 




16±f ± 2 


< 30 




16±6 


8.0-11.6 


K*+n° 


<31 






<31 


2.0-5.1 


K*°7t° 


<3.6 


<7 




<3.6 


1.8-4.4 



Table 3. Direct CP-asymmetry in B -> Kit, kk decays with = 40° ~ 90°, R b = \/p 2 + r] 2 = 0.38. Here we adopted = 1.3 GeV 
and ml = 1 .7 GeV for the PQCD results. 



Direct A C p(%) 


BELLE 


BABAR 


PQCD 


QCDF 


Charming Penguin ( \A cp ) 




77±27±8 


30±25±4 


16.0-30.0 


-6 ±12 


39 ±20 


n+n" 


-14 ±24^ 


-3±18±2 


0.0 


0.0 


0.0 


%+K- 


-7±6±1 


-10.2±5.0±1.6 


-12.9- -21.9 


5±9 


21 ±12 


% {) K- 


23 ± 11+ 1 


-9.0±9.0±1.0 


-10.0- -17.3 


7±9 


22 ±13 




7+y+i 

'-8-3 


-4.7 ±13.9 


-0.6- -1.5 


1±1 


0.0 



